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; MTHSHATANH returns the tke iy precision hyperbolic arctangent of the 
; H=floating precision argument. The call is standard call-by-reference. 
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-SBTTL DECLARATIONS 3; Declarative Part of Module 


: INCLUDE FILES: MTHJACKET.MAR 


: EXTERNAL SYMBOLS: 


; GLOBAL 
"SHOW BINARY CALLS g{QNDITIONALS, DEF INITIONS, EXPANSIONS 
*EXTRN MTHSK INVARG VARGH 


“EXTRN STNSeKD _ALOG_V 
-EXTRN MTHS$AB_H_FHT 


: EQUATED SYMBOLS: 


FWN (OOD NAMNE WN O OONAU EWR OOONAULE WhO 


0000 

B808 

Bn8 
8 00 

000 
8008 

000 
Bp08 

000 
0000 

000 

000 

000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 : 

0000 ; MACROS: 

0000 ; 

0000 

0000 MACRO OPDEF X, OP 

0000 OPDEF X “xO0@SH+0P,R'X,M'X 

0000 OPDEF ADDX3 *XO01Q@SH+OP,R'X,R'X,W'X 
0000 OPDEF xX “XO02@SH+0OP,R'X,M'X 
0000 OPDEF SUBX3 *XO0S@SH+OP,R'X,R'X,W'X 
0000 OPDEF LX *X04QSH+0P,R'°X,M'X 
0000 OPDEF MULX3 *XO5@SH+OP,R'X,R'X,W'X 
0000 OPDEF xX *X06@SH+0P,R'X,M'X 
0000 OPDEF DIVX3 “XO7@SH+OP,R'X,R'X,W'X 
0000 OPDEF CVTWXx *XOD@SH+O0P,RW,W'X 
0000 OPDEF POLYX *X15@SH+0P,R'X,RW,AB 
0000 OPDEF Vx “X7DFD,RH,WH ; MOVO 
0000 OPDEF MOVAX “X7EFD,AH,WL ; MOVAO 
0000 
0000 
0000 99 OPDEF 4H, <*X60FD> 
0000 -OPDEF ADDX, “XO 3608 8+*X60FD,RH,MH 
0000 -OPDEF ADDX3 *X0188+*X60FD,RH,RH,WH 
0000 OPDEF SUBX *X0208+*X60FD,RH,MH 

4 OPDEF SUBX3 ald a coutaee Ea’ aro 
000 OPDEF ULX *X0408+*X60FD,RH,MH 
0000 OPDEF MULX3 as Ce iy 9 oh Sa a 
0000 OPDEF DIVX *X%06a8+*X60FD,RH,MH 
0000 OPDEF DIVX3 *X07a8+*X60FD,RH,RH,WH 
0000 OPDEF CVTWX “X0D88+*X60FD,RW,WH 
0000 OPDEF POLYX “X15a8+*X60FD,RH,RW,AB 
0000 -OPDEF MOVX “X7DFD,RH,WH ; MOVO 
B00 -OPDEF MOVAX “X7EFD,AH,WL ; MOVAO 
000 100 
00000007 0000 £101 F_EXP = 7 ; Bit offset to exponent 
00000000 4 198 X_-EXP = 0 : Bit offset to exponent 


—> te 
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HATANH.MAR:1 99° (3). 
; PSECT DECLARATIONS: | 


03 

5 

8S ‘ -PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 
38 3 program section for math routines 
10 | 
11 
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; "FUNCTIONAL DESCRIPTION: 

HATANH = H=floating precision floating point function 

> HATANH(X) is computed using the following approximation technique: 
If iXi >= 1.0, error. Otherwise 

Let (14#X)/(1-x) = f * (2*n), where 1/2 <= f < 1 


If * Mi 4 tre than or equal to 1 then 
et 1 and Fl = 2ef. 
- 


set N =n and F = f. 


If iF = 1! < 2**-5 then 
2tatanh(X) = Neln(2) + W + WeP(W), 
mere t UCLSED/ (=P) )aoeen - 
and P is a polynomial of pm kek, F=5,D=9. 


OWONAO UE WN" OOONAUS WN OOon 


Else 
“obtain FHI (roughly squel to F) from 7 Loo 
2eatanh(X) = bn€(14+xX)/(1-x)) = Neln(2) + intend + 2*#Q(2*2), 
where Q is ? 9 poryosniat ef of ureere F=2,D=5, 


where Z = 
where F = (Deen #1930 712%) 


tet te thet tee te tee NNN NINO OO 


WR =O OONA UE WIN OOO NOUES WN $0 OONOU ES Wh 


Z is computed by: 
zZ = (Xx =D) /¢ =X*D) 
where = FHI *2eeN 

a = (¥=1)/(Y¥+1) f 

Note + Z may . — uted in a variety of ways 
Z = C(14#xX) = YJ/0(14xX) + ¥e(1-x)I 
Z=(1+X-e- Y ; Fn 1+X + Y¥ = X*¥) 
z=(il-Y¥*x + /L1 +Y¥+*+Xe- ah 
Z = ((1-¥) + XeC1ey) /UC14#¥) + Xe(1-¥) 


NOTE: The qyent tt tes Ln(A=FHI) and Ln2 are used in the above 
equet’ ons in two parts - a high part (containing the 
h hk ore bits) and a low pert (containing the low 
er bits. In the code the h r and low parts of the 
constents are indicated by a _H 
tively. The values were chosen such Ehat N*LN2_HI 
LN_FH “Hl is exactly representable. 


CALLING SEQUENCE: 
MTHSHATANH(hatanh.wh.r, x.rh.r) 
INPUT PARAMETERS: 


Y ; Address to store result 
x = 8 ; Contents of x is the argument 


; IMPLICIT INPUTS: none 


and _LO suffix coeges~ 
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oOating Precision Hyperb 6=-SEP 
OUTPUT PARAMETERS: 

VALUE: 
; IMPLICIT OUTPUTS: 
COMPLETION CODES: 
SIDE EFFECTS: 


none 
none 


; Associated message is: 
; CHFSL_MCH_RO/R1. 


preserves enables across the call. 


whether enabled or not. 


ERR: BRW ERROR 
ENTRY MTHSHATANH, ACMASK 
MOVX @X(AP), RO 
MOVX RO, R6 
SUBX3 R6, S*#1.0, -(SP) 
BLEQ RR 
ADOX $*#1.0, R6 
BLEQ RR 
CVTHF (SP)+, R4 


CVTHF R6, RO 


DIVF2 R4. R6 

MOVAB G“MTHSSAB_ALOG_V, R10 

ADDL2 (R10), R1 

BICWS «#laF EXP=1 R6, R5 

SUBW 36. #* 4000, R5 

BLEQ NEG_EXP 

SUBW Ss # 1aF_EXP, RS 

SUBW SRS, RG 

MOVZBL R6. R6 

IF NE F_EXP=X_EXP 

Divwe #13<F“EXP=X"EXP>, R5 

CVTWX RS, =(SP) 

CVTBL (R10)CR6], R10 

BLSS LN PLUS aw 

MOVAX G*AMTAS$SAB-H_FHICR10I, 
vx (R10)+, RB 
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H-floating precision hyperbolic arctangent of the argument 


Signals: MTHSK_INVARGMAT if iX! >= 1.0 with reserved operand in RO (copied to 
; the signal mechanism vector CHF$L_MCH_RO/R1 by LIBSSIGNAL). 

“Floating Overflow in math Library'’. Result is 
; reserved operand -0.0 unless a user supplied (or any) error handler changes 


NOTE: This procedure disables floating point underflow and integer 
; overflow, causes no floating overflow or other arithmetic traps, and 


; Note: This routine is written to avoid causing ony integer overflows, 
; floating overflows, or floating underflows or divid iti 


e by 0 conditions, 


standard call-by-reference entry 
gisetie DV (and FU), enable IV 
= arg 


(SP) = 1=-x 
ATANH(X) is not defined for X>=1 


ATANH(X) is not defined for X<=-1 


R6 = approximation to (1+#xX)/(1-x) 


R10 = address of ALOG table 
Biased exponent 

R5 = Unbiased exponent 

Branch to processing for n=<0 
Exponent is positive, RS =N=n- 1 
R6 = F = 2f 

R6 = index into ALOG table 


Shift RS to scale X-floating 


Push N onto the stack 

R10 = offset into FHI tables 

Branch to handle F close to 1 
; R10 = Address of FHI 

R6 = FHI 


1-0 
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MTHSHATANH ;_H=floating Precision Hyperbolic Arctan 16-SEP-1984 01:34:37 VAX/VMS Macro v04-00 Page 
2-002 THSHATANH © H- floating io na be. Hyperb aes 913g: 27 MTHRTL.SRCJMTHHATANH.MAR; 1 . &) 
74 a. 
74 § ; Compute Z = (F = FHI)/(F + FH 1} 
74 : 2 = b(t sx) = Ye(4 mary (14x) + Y¥#(1-x) 
74 4; Z=(1+X- Y¥ + XxevyJ/[1 +X + Y = XeY 
i ; 3 where Y = FHI*2**N, roughly equal to (1+#X)/(1-x) 
56 55 AO 0274 287° ADDW. RS, -R6 ; R6 = FHI * geen = SFHI 
7E 08 56 63FD 77 88 SUBX3 R6, $*#1.0, -(SP) ; ($P) = 1 - 
56 08 60FD Oe7c 289 ADDX $°#1.0, R6 : R6 = 1 + SFH ih 
6E 56 66FD 0280 290 DIVX. R6, (SP) + (SP) = (1=SFHI)/(1+SFHI) = D 
56 6€ 50 61FD 84 91 ADDX3 RO, (SP), R6 ; R6 =D +X 
6E 30 64FD 89 3¢ MULX RO (SP) : (SP) =D * X 
6— 08 60FD dD 29 ADDX $441.0, (SP) : (SP) = 1 + Dex 
56 BE 66FD 0 4 36 DIVX (SP)+, R6 > R6 = (D4X)/(14D*xX) = 2 
8 32 38 ; Compute Z2**2, P(Z**2) and Z2*P(2**2) 
50 56 56 65FD 0295 98 MULX3 R6, R6, RO ; RO = Z**2 
FEBF CF OA 50 75FD 029A 99 POLYX RO, #LOGLEN2, LOGTAB2 ; RO = P(Z**2) 
50 56 64FD OgAl 99 MULX R6, RO 3 RO = Z2*P(Z**#2) 
BgA2 8 ; Compute B = N*LN2_LO + LN_FHI_LO + Z*P(Z*Z) 
56 FD65 CF 6€ 65FD 02a5 04 MULX3 (SP), LN2_LO, R6 ; R6 = N*LN2_LO 
56 8A 60FD O2AC 305 ADDX (R10S+, Rb > R6 = N®LNe- to + LN_FHI_LO 
50 56 60FD 9280 306 ADDX R6, RO ; RO =8 
ase? 38 ; Compute A = N*LN2_HI + LN_FHI_HI and ALOG(X) 
56 FD46 CF 8E 65FD 0284 10 MULX3 (SP)+, LN2 LHI, R6 : R6 = N*LN2_HI 
56 6A 60FD 0288 11 ADDX (R 10) R6 ; RE= A= N*LN2_HI + LN_FHI_HI 
50 56 60FD O2BF 31¢ ADDX RG ; RO = A+ B= ACOG(X) 
50 01 A2 6 (3) 31 SUBW2 RiaKe EXP, RO ; Divide by 2 
04 BC 50 7DFD C6 14 MOVX RO, SYC(AP ; Store result 
04 02CB 15 RET 
02cC 16 
02cC 17 NEG_EXP: 
5 55 2 O2cc 18 SUBW SS séRS,, -R6 ; R6 = F = 2f 
56 56 QA O2CF 319 MOVZBL R6, R6 : R6 = index into ALOG table 
00000007 0 D 30 IF NE, F_EXP=X_EXP 
55 0080 8F Ab p 1 Divwe #1a<FTEXP=X"EXP>, RS s-;_: Shift RS to scale ies. 
7E 55 6DFD § D7 § CVTWX RS, -(SP) ; Push N onto the st 
SA 6A46 O98 9 DB 4 VIBL (R10)CR6], R10 : R10 = offset into rit tables 
DF 5 LN_1_PLUS W: 
65 19 DF § eCss PLUS Branch to handle F close to 1 
5A 00000000'GF4A 7EFD E1 MOVAX GATRSSAB_H.FHICRTON, nib ; R10 = Address of FHI 
56 6A 7DFD EA § MOVX 10), R6 R6 = FHI 
EE 0 ; Compute Z = (F = FHI)/(F + FHI) 
EE 13 Z = C(14x) = Ye eae: . votts =a] 
Ee ¢ : 2=f1 +X -Y¥ + xeV)/ 
43 ? : where Y = FHI*2**N, roughly equal to C070 
56 55 Ad EE 5 ADDW = RS, -R6 : R6 = FHI * geen = SFHI 
7E 08 56 63FD F $ SUBX3 R6, $*#1.0, =(SP) : (SP) = 1 = SFHI 
6 08 60FD 02F6 ADDX  $*#1.0, R6 > R6=1 + SFAL 
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2-002 MTHSHATANH : H-floating oP ec sion Hyperb o78Fb= 1986 913d ib6 EMTHRTL. SRCIMTHHATANH. MAR: 1 _ 3) 
bE 66FD O2FA DIVX. R6, (SP) : (SP) = (1-SFHID/(1+SFHI) = 
56 COG 61FD OoF 38 ADDX3 nS. (SP), R6 : R6 =D +X : sircse te 
&E 64FD 0 40 MULX : (SP) =D * X 
6€ 08 60FD 4] ADDX  $*#1.0, (SP) : (SP) = 1 + Dex 
56 BE <66FD 0B 4g DIVX (SP)+, R6 + RO = (D+X)/(14D#x) = 2 . 
. te ; Compute Z2**2, P(Z**2) and Z*P(Z**2) | ° 
50 56 56 65FD O30F 346° MULX3 6, R6, RO 7 RO = 2*22 1 
FEGS CF =0OASs«SS Tas 0314 23 POLYX RO. #LOGLEN2, LOGTAB2 ; RO = P(Z*#2) 2 
50 56 64FD 0 18 48 MULX R6.R : RO = Z2*P(Z**2) - 
aM 29 : Compute B = N*LN2_LO + LN_FHI_LO + Z*P(7Z*Z) LON 
56 FCEB CF 6€ 65FD O31F 52 ° MULX3 (SP), LN2_LO, R6 : R6 = N*LN2_LO 
56 7A 60FD 0326 2g ADDX <= (R16), RB ; R6 = N*LN2-LO + LN_FHI_LO 
50 56 60FD 03 A 34 ADDX 6, RO > RO=8 
03 : 36 : Compute A = N*LN2_HI + LN_FHI_HI and ALOG(X) 
56 FCCC CF 8 6BE 65FD «(03 58° MULX3 (SP)+, LN2_HI, R6 > R6 = N*LN2_HI 
56 7A 62FD «(03 i 59 SUBX =(R105, R6” ; R6 = A = NELN2_HI + LN_FHI_HI 
50 56 60FD 0339 360 ADDX R6, R >: RO = A+B = ALOG(X) 
50 01 A2 0330 361 SUBW2 #1aX_EXP, RO : Divide by 2 
04 BC 50 7DFD bae0 385 MOVX RO, SYC(AP) ; Store result 
04 034 6 RET 
b3ce 308 
0346 66 : Special logic for F close to 1 
gee 
0346 369 LN_1_PLUS: 
56 08 50 63FD 0346 70 SUBX3 RO, S*#1.0, R6 : R6 = 1-x 
33 f re 4 ai ah 3 Determine which way to calculate Ww 
56 10 56 67FD baa 8 DIVX3 R6, $*#2.0, R6 > R6 = 2/(1-Xx) 
36 08 62FD 0394 i SuBx 5481.0, R6 ; R6 : (07-0 
P cate 
56 08 62FD si 26 SuBx 4#1.0, R6 > R6 = W 
56 «550 98 67FD 8 61 28 10S: DIVX3 R6, RO, R6 > R6 =X / (1-x) 
36 gi Ad 366 7 ADDW #1aXx_EXP, R6 + RO = W = 2eX/(1=X) = (14X)/ (12K) = 1 
FCBO CF 13 28 Tarp $9 30 20$: pou R6, PLOGLEN1, LOGTAB1 ; Re = Vo ae 
7E FC96 CF 6E 65FD Ni : MULX3 (SP), LN2_LO, -(SP) : (SP) = N*LN2_LO 
50 8 60FD 0378 ADDX  (SP)+, RO > RO = Nstg tO + Wea(W) 
50 56 60FD 037F 4 ADDX R6, R : RO = N*LN2-LO + WeQ(W) + W 
6E C78 CF 64FD x 5 MULX LNO_HI, (SP) : (SP) = N*LN2_HI 
30 BE 60FD 9 : ADDX  (SPT+, RO : RO = ALOG(X) 
0 gi A2 SUBW2 #lax_EXP, RO : Divide by 2 
04 BC 0 7DFD p 8 MOVX RO, SY (AP) ; Store result 
396 (390 ; ™ 
96 391 ; X <= 0.0, signal error | 
9A 3 38 ERROR: MOVZBL #MTHSK INVARGMAT, =-(SP) ; condition value 
79 9A 94 ASHQ #15, #T, RO ; RO = result = reserved operand -0.0 
} 
| 
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MTHSHATANH 
2-002 MTHSHATANH = H-floating Precision Hyperb 6-SEP=1984 11:24 


52 7C «(039 95 CLRQ ss R2 
A 96 : 
ve : 
00000000'GF 01 FB AQ 399 CALLS #1, G*MTHSS$SIGNAL i § 
04 BC 50 7DFD 03A7 431 MOVO RO, @Y(AP) : 
04 O3AC 4 é RET : 
40 
404 END 


4:37 VAX/VMS Macro V04-00 Page 
746 CMTHRTL.SRCJMTHHATANH.MAR; 1 


goes to signal mechanism vector 
CHFSL_MCH_RO/R1) so error handler 
_‘- modify the result. 
“oe error and use real user's PC 
ndependent of CALL vs JSB 
ee result 
return 


10 
(5) 


————— 
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MTHSHATANH oat eh 1: 4:37 VAX/V 
Symbol table 6-SEP=-1984 11:24:46 (CMTHR 
ACMASK = Q00007FC 
ERR 00 Sei? & 4 
ERROR 4 88 B 1 
F_EXP = 0000000 
LN2_HI sid os st R 01 
LN2_LO 464 10 R 01 
LN_T_PLUS 0000346 R 8 
LN_1 PLUS_W si it DF R 1 
OGLEN1 = 00000013 
LOGLEN2 = QOQO0000A 
TABI 00000020 R 01 
LOGTAB2 00000160 R 01 
MTHSSAB_ALOG V  teeeneee =X 00 
MTHS$AB_H_FHT eererere = X 00 
MTHSSSIGNAL eeneeeee = =X 00 
MTHSHATANH 00000213 RG 01 
MTHSK_INVARGMAT eeeeneee x 00 
NEG_EXP QOO0002CC R 01 
= 00000008 
X_EXP = 00000000 
Y = 00000004 
$e sr en wares wn on oe + 
! Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
. «ABS. 00000000 < 0.» 00 ¢ OO.) NOPIC USR CON ABS LCL NOSHR NO 
_MTHSCODE QOOOO3AD (¢ 941.) O1¢ 1.) PIC USR CON REL LCL SHR 
pow wr ee wwe rer anw ew ewe nero 
! Performance indicators ! 
teow ce eee eee ee eee een eee eee 
Phase Page faults CPU Time Elapsed Time 
Initialization — 34 00:00:00.1) 00:00:01.13 
Command processing 134 00:00:00.66 00:00:04.89 
Pass 101 00:00:01.55 a BE rH 
Symbol table sort 0 00:00: 8-9 Ba BR: 1 
Pass 2 % Bp: 0:00.99 0:00:03.50 
Symbol table output :00:00.0 0:00:00.0 
Psect synopsis output 3 :00:00. 8:08: 8 0 
Cross-reference output :00:00. :00:00.0 
Assembler run totals 366 00:03. 0:00:14.88 


The working set Limit was 900 pages. 

7868 bytes (16 pages) of virtual memory were used to buffer the intermediate gefe. 

There were 10 pages of symbol table space oh ecatee to hold 21 non-local and 2 local symbols. 
464 source lines were read in Pass 1, producing 15 object records in Pass 2 

3 pages of virtual memory were used to define 2 macros. 


Page 


EXE NORD NOWRT NOVEC BYTE 


RD NOWRT NOVEC LONG 


VAX=11 Macro Run Statistics 


Macro Library name Macros defined 
“$255$DUA28: CSYSLIBISTARLET.MLB;2 0 
0 GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 
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MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHATANH/OBJ=OBJ$:MTHHATANH MSRC$:MTHJACKET/UPDATE=(ENHS: MTHJACKET) MS 
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